Quintom cosmologies with arbitrary potentials 
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In this paper we investigate quintom cosmologies with arbitrary potentials from the dynamical 
systems perspective. The dynamical systems analysis is complete in the sense that it includes the 
asymptotic regime where both scalar fields diverge, which proves to be particularly relevant in 
connection with the existence of tracking phases. The results of the present study indicate that the 
existence of phantom attractors is not generic: for quintom models there may exist either de Sitter 
attractors associated with the saddle points of the potential, or tracking attractors in the asymptotic 
regime where the scalar fields diverge. 
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I. INTRODUCTION 



Observational cosmology supports the claim that our 
universe is currently undergoing accelerated expansion 
[3; S 9 • Within the framework of Einstein's theory of 
general relativity, the arguably only possible explanation 
to this phenomenon relies on the existence of an enig- 
matic component of the universe acting as a repulsive 
force, the so called dark energy (DE). 

Dark energy is one of the most debated topics in the 
era of precision cosmology (see [H, H[ for recent reviews). 
A cosmological constant (CC) is the simplest candidate 
to account for the DE, but this proposal suffers from the 
well-known cosmological constant problem according to 
which there exists a (yet) unexplained extraordinary dis- 
crepancy (of about 120 orders of magnitude) between the 
(cosmologically) observed value of the CC and the value 
predicted by quantum field theory. For this reason, alter- 
native routes has been proposed support (see again 
quintessence, phantom dark energy, tachyon fields, Chap- 
lygin gas, and recently, quintom dark energy. The latter 
alternative has aroused interest due to a debate about 
the possibility that the DE equation of state (EOS) pa- 
rameter w(z) crosses the phantom divide line w = — 1 
at recent redshifts. 1 If this possibility were eventually 
corroborated, then one could think of two possible ex- 
planations/implications: either the DE consists of multi- 
ple components with (at least) a non-canonical phantom 
component, or general relativity has to be modified on 
cosmological scales Q ■ The possibility to deal with sev- 
eral scalar fields is supported, besides, by the effective 
(low-energy) theory of the bosonic string where several 
scalar fields like the dilaton, the axion and other moduli 



fields can be found [8(- 

If dark energy becomes phantom-like (phantom dark 
energy has an equation of state (EOS) parameter w 
smaller that —1) at recent times then, typically, our uni- 
verse (if ever expanding) can evolve to a catastrophic 
singularity in the near future, characterized by diver- 
gences of the scale factor, the Hubble expansion scalar 
and its time derivative. It is the case if the dominant 
fractional energy density is phantom like. The phantom 
dominated late time solution appears in the context of 
quintom cosmologies with exponential potentials when 
one resorts to standard dynamical systems techniques as 
it has been shown in 0, [l^] and . The quintom is a 
hybrid of a quintessence component, usually modeled by 
a real scalar field that is minimally coupled to gravity, 
and a phantom field: another real scalar field -minimally 
coupled to gravity- with negative kinetic energy [9|, 1 10 . 

The aim of the present paper is, precisely, to provide 
further elements to support the interest on quintom DE 
models. We do so by studying from the dynamical sys- 
tems perspective quintom cosmologies with arbitrary po- 
tentials (other than exponential ones). In this regard, 
our study is complementary, in some sense, to that of 
references [§, [l(J El ■ As customary, this kind of analy- 
sis relies on an appropriate choice of phase-space coordi- 
nates, which builts partially on those made in precedent 
references [TTI. [27j. 

Perhaps, the most important conclusion of the present 
research is confirming the result of reference fllj ] (the ex- 
istence of phantom attractors is not generic) in the more 
general setup of quintom DE with arbitrary potentials 
(other than exponential ones). 3 Peculiarly, in our frame- 
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1 The quantum stability of a w < — 1 phase of cosmic acceleration 
was analyzed in [f|. 



2 Quintom-like behavior (with w = — 1 crossing) has been found 
also in the context of holographic dark energy 24, 25, 26]. 

3 The same conclusion was reached in reference ll-ll . where a 
hessence model with exponential and inverse power law poten- 
tials were studied from the dynamical systems view point. 
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work, de Sitter (w = —1) attractors can be associated 
only with saddle points of (the natural logarithm of) the 
potential. 

The paper has been organized as follows. In section 
II we give the details of the dark energy model under 
investigation. In section III we turn to perform the dy- 
namical study of this quite general cosmological model 
for finite values of the component scalar fields. For illus- 
tration purposes, we choose a class of potentials which 
was derived in the literature on inflation. Then, and in 
order to complement the dynamical system analysis car- 
ried out in section III, we make in section IV a different 
choice of phase space variables relating the scalar fields, 
which allows to study the asymptotic regime where both 
component scalar fields diverge. Conclusions are given in 
section V. Finally, to make the paper self-contained, we 
have added an Appendix with some recipes on dynamical 
systems and its applications. Throughout this paper we 
use natural units. 



II. THE MODEL 

We will investigate the evolution of a spatially flat 
Friedmann-Robertson- Walker (FRW) spacetime fuelled 
by dark matter (dust) with energy density p m , and by 
quintom dark energy with energy density p de and pres- 
sure pde given respectively by 

Pde = \^> 2 - ^p 2 + V(cf>,tp), p de = ujp de , (1) 

where for our quintom model, the scalar field tp plays the 
role of a phantom field. In the above equations 

^-p 2 -2V{^p) 

w = - (2) 

<f) 2 -p 2 + 2V(<f>,p) 

is the equation of state (EOS) parameter of the dark en- 
ergy. Thus, we are assuming our dark energy comprises 
the contribution of two scalar fields <p and ip interact- 
ing through their potential energy V((j>, tp), and satisfying 
separate energy conservation equations, namely: 

4> + 3Hj, + d^V{(p, tp) = 0, (3) 
tp + 3Htp - d^V^, tp) = 0. (4) 

We adhere here to Einstein's theory of general rela- 
tivity so, as often done, dark energy will be assumed to 
be minimally coupled to matter. With our assumptions, 
the cosmological equations are the standard Friedmann 
equation 

3H 2 = p. m + p de , (5) 
the continuity equation 

p rn + 3Hp m = , (6) 



and the conservation equations Eqs. ([3l 0J. Here and 
throughout overdots denote differentiation with respect 
to cosmic time t, H = a/a is the Hubble factor, and a is 
the scale factor. 

Combining Eqs. $S§§ one obtains the evolution equa- 
tion for H: 

- 2H = p rn + p de + p de . (7) 

Equations ([3][7]) are the basic set required for formulat- 
ing a dynamical systems investigation of the evolution of 
these cosmological models. 



III. PHASE-SPACE 

It is well-known that information about the evolu- 
tion of cosmological models can be retrieved using dy- 
namical systems as a tool. Particularly, the asymptotic 
(and sometimes the intermediate) behavior of cosmolog- 
ical models are closely related to concepts like past and 
future attractors (28l. |29||. Our purpose, thus, is to take 
advantage of that fact, and apply it to the model under 
consideration. 

We want to construct a dynamical system with the 
following properties: 

1. It is autonomous. 

2. It has the lowest possible dimensionality. 

3. It describes a physically interesting and reasonable 
configuration within the framework of quintom cos- 
mologies. 

In order to construct a dynamical system with these 
properties one of the first steps is to introduce a set of 
convenient (in most occasions expansion normalized (29[) 
variables which allow rewriting the conservation equa- 
tions and the evolution equation of H as a dynamical sys- 
tem subject to a constraint arising from the Friedmann 
equation ([5]). Since we want to study arbitrary poten- 
tials (as arbitrary as possible, but with the exclusion of 
exponential like potentials), a minimum of six variables 
is needed to construct an autonomous dynamical system. 
However, the problem gets somewhat simplified -at least 
concerning its dimensionality- if we use the constraint 
arising from Eq. |5]) as a definition for one of the vari- 
ables. Hence, the minimal number of dimensions of our 
model is lowered down to five. 

The next step in the study of the evolution of our dark 
energy model viewed as a dynamical system is to find its 
fixed (or critical) points. The stability of the fixed points 
is then analyzed by studying the linearized dynamical 
system obtained by expandin g th e evolution equations 
about those fixed points (see [28l [3(| for seminal refer- 
ences and the Appendix for useful recipes). 

Turning back to the first step, it follows that one 
can present Eqs. ([SK]) in the form of a dynamical 
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system by making the following choice of variables: 

(x^, x v , y,z, <j>,ip), where 



X<f> 



V6H 7 



V3H' V3H' 



(8) 
(9) 



and <f), ip are the scalar field variables. This choice of 
phase-space variables renders the Friedmann equation as 



2,2,2 

x v + y +z 



1. 



(10) 



The constraint ([10]) follows from Eq. (J5J) and it allows, 
as we have noticed before, to consider only the evolution 
of those five variables other than z (the evolution of which 
is being determined by the evolution of the former ones). 

In what follows we will restrict ourselves to y > and 
to H > since the main concern of this paper are ex- 
panding universes. Combining expressions ([3][9]), the fol- 
lowing evolution equations are obtained: 



V6 



y 2 d (jl \nV + {q-2)x (j> , 



= U^y 2 d v> lnV+(q-2)c 



(11) 



(12) 



— 1 + q 

3 V H 2 

V6 



3 

V6 



lnV + a ¥ v lnV))y,(13) 
(14) 
(15) 



Here primes denote differentiation with respect to a new 
time variable r = log a 3 , and q = —da/ a 2 stands for the 
deceleration factor. Explicitly, 



In general, the phase portrait can be split into two re- 
gions depending on the values of </> and ip. One region 
is defined by the condition that both scalar fields are fi- 
nite. This region comprises matter dominated and de 
Sitter solutions. In the second region both scalar fields 
diverge. As we shall see, it is the latter region the one 
that may accommodate phantom solutions whereas the 
former cannot, that, is, phantom behavior only appears 
in the regime where the scalar fields blow up. In this 
regard we have to notice that the above choice of phase- 
space variables (I1HI15[) is not adequate to study the re- 
gion where scalar fields diverge, so a different choice of 
variables is necessary. This case (infinite (4>,tp)) will be 
studied separately in section IV. 



Matter dominated solutions 



The dynamical system (|1HI15P admits the biparamet- 
ric class of non- hyperbolic critical points (x$, x v ,y) = 
(0,0,0), for arbitrary (<j>,<p) G R 2 . These critical points, 
which we will denote generically as O, represent matter 
dominated cosmological solutions. Typically, the scale 
factor and the matter energy density evolve as a a t 2 ^ 3 
and p m oc t~ 2 respectively. This class of solutions is rel- 
evant from the end of the epoch of radiation domination 
till the very recent epoch when the dark energy contri- 
bution equated the contribution from the dark matter. 

Making a local linear analysis we obtain that the eigen- 
values of the linearization around these critical points 
(with (j) and ip fixed) are (-1/2,-1/2,1/2,0,0). Due to 
the existence of two null eigenvalues, as already said, 
these critical points are non-hyperbolic so that the linear 
analysis is not conclusive in this case. Alternatively, one 
can determine analytically the stable, unstable and cen- 
ter manifolds of those critical points (see the Appendices 
A and B). 



«-i(3 



(2 2 2 



1) 



(16) 



The evolution equations of variables x^,, x v , y, (f> and (p 
form a 5D dynamical system defined on the phase space 



* = {(a?*, x v , y) : < x\ - x% + y 2 < 1} x 

{(^)GR 2 }. (17) 

The phase space *f> is unbounded because the variables 
x v , (j), tp can reach infinite values. However the projec- 
tion of the phase space into the subspace {x r j >1 x v , y) is 
limited by the hypersurface defining the hyperboloid. 



B. de Sitter solutions 



The dynamical system (|1HI15P admits solutions where 
(x^,x v ,y) = (0,0,1), for each (</>*, ip*) 6 R 2 such that 
d^\nV{<j)*, V *) = d v lnV(^,p*) = 0, that is, (</>*, tp*) 
are the stationary points (extrema or saddle points) of 
InV. 5 

For fixed <fi* and ip* the corresponding solution is dom- 
inated by the potential energy of the quintom field, which 
means it is a de Sitter solution, so we will denote this class 
of fixed points as dS. The scale factor evolves according to 

the following exponential law a oc exp y/V(cj>*, ip*)/3t 

These solutions are expected to be important at late 



4 Notice that the equations I|1H15|I are invariant under the change 

y -* -V- 



5 The stationary points of InV coincide in general with the sta- 
tionary points of V provided these are not zeros of the potential. 
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times when the potential energy mimics an effective cos- 
mological constant. 

The eigenvalues of the linearization around the corre- 
sponding fixed points are 



Ai = -1, 




(18) 
(19) 

(20) 



where 



Ai = 



d 2 \nV d 2 \nV 

dip 2 dcj) 2 



d 2 \nV d 2 \nV 

d<j> 2 dip 2 

<9 2 lnV x 2 



(21) 



6* ,ip=ip* 



dipdcj) 



(22) 



If the point (</>*, ip*) is an extremum of InV then the 
quantity A defined by 



A = 



d 2 \nVd 2 \nV fd 2 \nV 



dip 2 dcj) 2 



dipdcj) 



(23) 



,<p=v>* 



is positive; if it is a minimum (maximum) then 
d 2 \nV/d(t) 2 > (< 0) and d 2 \nV/dip 2 > (< 0). If 
A < 0, then {(f)* , p*) is a saddle point of In V . 

From the former analysis it follows that the dynami- 
cal behavior of the critical points in 'J depends on the 
classification of the critical points of In V: ((f)* , ip*). 



1. Dynamical behavior of the critical points in ^ associated 
with the de Sitter phase 

If ((f>*, ip*) is an extremum of \nV(<f) 1 ip), then the con- 
dition A > implies A2 > A 2 . In this case the corre- 
sponding critical point in is a saddle; i.e., all the real 
parts of the eigenvalues of the matrix A evaluated at the 
corresponding critical point are different from zero and at 
least two of them have different signs (see the Appendices 
A and B). 

If ((f)* , ip*) is a saddle point of In V((f>, ip), then the con- 
dition A < implies A2 < A 2 . In this case the corre- 
sponding critical point in the phase space \1/ is an attrac- 
tor in the following three cases 6 : 



• Case i) All the eigenvalues are negative reals. This 
is the case if 

< A 2 < 9/16 and - 3/2 + yfK~ 2 < Ai < -^/A^. 



• Case ii) A* are complex conjugate eigenvalues and 
Xf are negative real eigenvalues. This is the case if 



< A 2 < 9/16 and - VA? < 3/2 + Aj < VA2, 



or 



A 2 > 9/16 and - 3/2 - VA^ < Ai < -A/A2. 

• Case iii) A^ and A 3 are respectively complex con- 
jugate eigenvalues. This is the case if 

A 2 > and Ai < -3/2 - y^. 



C. Example: hybrid inflation potential 

To illustrate the former result we will consider a two 
scalar fields model with potential 



V(*, V ) = \m*f + lgW + ~ A A ^ )2 , (24) 

where m, g, M, A are real-valued non-negative con- 
stants. 7 

The potentia l (1241) has been formerly considered in 
USES HUM [H[ in the context of hybrid in- 
flation. Since our context is late universe, the virtues of 
this potential regarding (early time) inflation are not of 
interest to us, we have selected it because its functional 
form is suited to illustrate some of the peculiar features 
of the late time evolution of our quintom dark energy 
models. 

It is straightforward to show that In V((f>, ip) has 
a unique real-valued stationary point at (<f)*,ip*) = 
(0,0) which is a saddle point of InV since A(0,0) = 
-16A 2 m 2 /M 6 < (see eq. This means that the 

de Sitter phase is an attractor in ^. Actually, since the 
function hiV with V given by ([24]) has only a saddle 
point (it has no minimum at all), for a quintom model 
the critical point (in the phase space \1/) associated with 
(((>*, p*) is stable in the following three cases as follows 
from the former analysis (see cases i), ii), and iii) in the 
former subsection): 



The critical point is an the attractor if all the real parts of the 
eigenvalues of the matrix A evaluated at the corresponding crit- 
ical point are negative (see the Appendices A and B). 



7 Quintom cosmologies with polynomial potentials were investi- 
gated in [31I , |33 | , where only dark energy (no CDM) was consid- 
ered. 
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(a) 




(b) 




(c) 



FIG. 1: Projections of the phase-space trajectories for dif- 
ferent values of M, m, g, and A: M = 4, m — 3, g = 2, 
A = 2 (a); M = 2, m = 1.5, 3 = 2, A = 1.3 (b) and 
M = 4, m = 5, g = 2, and A = 3.5 (c). The point dS rep- 
resents the de Sitter attractor. All the orbits portrayed have 
been truncated so that y > and they have been obtained by 
fixing the initial values for the scalar fields to cj>(r = 0) = 0.20 
and ¥>(t = 0) = 0.19. 



• Case i) This will the case (see Fig. QJ,) if either 

< M < to and < A < 3Af 4 /16m 2 
or if 

< to < M and < A < 3Af 2 /16. 

• Case ii) This case will occur (see Fig. |T|d) if either 
< M < m and 3M 4 /16to 2 < A < 3M 2 /16 

or if 

< m < M and 3M 2 /16 < A < 3A/ 4 /16to 2 . 

• Case iii) Finally, one will have this case (see Fig. 
Ht) if either 

< M < to and A > 3Af 2 /16 

or if 

< to < M and A > 3M 4 /16to 2 . 

This completes the study of the region with no diver- 
gent scalar fields. 



IV. EQUILIBRIUM POINTS AT INFINITE (0, tp) 

One of the most widespread view points within the 
framework of (DE) scalar field cosmologies is that only 
exponential potentials may be responsible for a tracking 
phase. For instance, exponential potentials have been 
considered as a source of interaction between scalar fields 
in conventional [13, EH S3, 511 and unconventional 0, 
Hoi . [TlT | cosmologies. However, in reference [l4|, it was 
shown that in the context of quintom-likc DE cosmologies 
(hessence models) with exponential and (inverse) power 
law potentials, some stable attractors can exist which are 
either scaling solutions or hessence-dominated solutions 
with EOS larger than or equal to -1. Similar results will 
be confirmed in what follows. 

In [ll[ a step forward was done with respect to [§, [l(| : 
The potential used in the former reference 



(25) 



leads to the existence of an scaling regime; i.e., it was 
proved that quintom cosmologies with exponential poten- 
tials were not devoid of asymptotic scaling phases. In 
the context of this study, the question that comes to 
mind is the following: can there exist tracking solutions 
in the quintom scenario with arbitrary potentials? This 
question is partially motivated by the fact that track- 
ing behavior may exist in single scalar fields models with 
potentials that asymptotically tend to exponential ones 
|27| . so it might well be the case that in our context 
a broad family of potentials satisfying some particular 
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requirement regarding their asymptotic functional form 
could display tracking behavior. 

The existence of tracking phases in the region of finite 
<t> and ip is precluded by our previous analysis. So one 
has to turn attention to the region where <f> = ±00 and 
<p = ±00. Thus, our purpose is to try and investigate 
the possible existence of scaling phases in quintom cos- 
mologies with arbitrary (other than the exponential) po- 
tentials, by concentrating on the region where the scalars 
fields diverge. This is to be done following the method 
developed in j22|. It must be stressed that our study of 
this region will be done in full detail, and even though 
emphasis will be given to tracking solutions, other possi- 
ble asymptotic states will be reported in case they exist. 
This analysis will be complementary to the dynamical 
study of section III. 

We start by making a different choice of phase space 
variables relating the scalar fields: u — v = ip^ 1 . 

The latter variables are complemented by x^,, x v , and y 
defined in Eqs. ©■ Making use of the required coordi- 
nate transformations the original system (|1HI15[) can be 
expressed as 



i ((-2 + q)x 4> + 3Ny 2 ) , 



3 ( (-2 + q)x v - 3My< 



-y\l + q-3[Nx4,+Mx v 
\/6 2 

r-W X 



V6 



V 2 X a 



(26) 
(27) 
(28) 
(29) 
(30) 



where V , N and M are arbitrary functions of u and v 
defined, respectively by: 



V(u,v) = V(u -1 ,u _1 ) 
\/6 (N, M) 



nd\nV odlnV 
u — - — ,v 



On 



(31) 
(32) 



We assume V is such that 



lim 

(u,u)^(0,0) 



2 dlnV 2 dlnV 
u — - — , v 



du 



dv 



V6(n, fh). (33) 



The requirement (|3"3"|) implies that the only essential 
feature of the potential is that it should behave asymp- 
totically as an exponential potential. This (very) general 
condition on the potential V, leads us to study the phase- 
space trajectories in the invariant subset (u, v) = (0, 0) 
by taking the limit (u, v) — > (0, 0) in the equations (|261 
I28p . By continuity, the reduced dynamical system in that 



invariant subset is 

4 = I ( 3 ™/ 2 + (l - 2)a*) , (34) 
x' v = ~(3hy 2 -(q-2)x v ), (35) 

y' = g(l + q - 3(ma;0 + nx v ))y (36) 

The equations (|34M36p are the same as in [ll|, but the 
subtle point here is that the dynamical system (|34H36j) 
describes the dynamical behavior of a quintom cosmology 
with an arbitrary (asymptotically exponential) potential 
in an invariant subset of R 5 . So the dynamics is richer 
than in Nevertheless, the results reported in the 

latter reference apply to our more general case. These 
are immediately summarized. 

As we have noticed before, the first approach to ana- 
lyzing the evolutionary behavior of a cosmological model 
from dynamical systems point of view is to find the crit- 
ical points of the system arising when normalized vari- 
ables are taken into account. Particularly at the invariant 
subset (u, v) = (0, 0) there may exist representatives of 
the critical points O, C±, P and T studied in [H|]. 8 ^ n 
Table [J the location, the existence and deceleration fac- 
tor of these critical points for fh > 0, n > and y > 
are presented. 

The existence of a scaling regime of is represented by 
the critical point T for values of m and n such that 
m 2 — n 2 > 1/2. Whenever it exists, it represents a attrac- 
tor solution in which quintom dark energy tracks matter 
(the equation of state of the quintom fluid is dust-like; 
i. e., w = 0). This late-time asymptotic state rep- 
resents a decelerating cosmological model in which the 
fractional energy densities of matter and dark energy 
are proportional. Particularly, f^ e and Q m are given by 
fl de = (2m 2 -2n 2 )- 1 , and tt m /fl de = -l+2m 2 -2n 2 > 0. 
The dynamical character of T depends on the value of 
fide- it is either an stable focus or an stable node if 
Qde < 7/8, or 7/8 < VLde < 1, respectively. 

The fixed point P represents a solution in which quin- 
tom dark energy dominates over matter (the equation of 
state of the quintom fluid corresponds to a fluid that red- 
shifts faster than dust). When such point exists it is an 
attractor only when its existence precludes that of the 
point T. This late-time asymptotic state does not neces- 
sarily represents an accelerating cosmological model, that 
depends on the quantity fh 2 — n 2 . Given its eigenvalues 
structure (see table H]), this solution is either a saddle or 
a stable node. The accelerated solutions associated with 
this fixed point can provide a good representation of the 
presently observed universe. 

It is noticeable that in the regime (0 ^ 00, ip — > 00) 
phantom solutions may exist as well. Actually, a po- 
tential satisfying (f33|) allows P to be characterized by 



8 For fixed points of class O the interpretation in Section IIIA is 
applicable as well. 
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TABLE I: Location, existence and deceleration factor of the critical points for m > 0, h > and y > 0, in the infinite (<f>, <p) 
region. 



Name x,p x v y Existence 



O All rh and n 



C± iyl + i* 2 x* All m and n 2 

P rh -h Vl - m 2 + n 2 ra 2 - n 2 < 1 -l + 3(ra 2 -n 2 ) 

~ m n 1 -2-2^-, /0 1 

m — n > 1/2 



2(m 2 -n 2 ) 2(m 2 -n 2 ) 27m 2 - n 2 ' 2 



w < — 1 (i.e. phantom dark energy) if rh 2 < n 2 . Other 
possibilities include the cases w = —1 (i.e. de Sitter vac- 
uum dark energy) if rh = n, and —1 < w < —1/3 (i.e. 
quintessence dark energy) if < rh 2 — h 2 < 1/3. These 
are the solutions which the interest on quintom cosmolo- 
gies is grounded on. 

In general the potential V responsible for such a scaling 
regime can be reconstructed from functions N and M 
satisfying: 



v~ 2 —— =u~ 2 — , hm (N,M) = (n,m) 

OU OV (u,ti)->(0,0) 



as 



V [u, v) = Vo exp J a/x x 



( V6M(u, v) f u V6d v N{n, v) 

P / : 1 ^ ./, r 



dn dv. (37) 



Note that 



d 2 \nV V6dM 



dudv 



du 



Finally, this region of the space phase admits non- 
isolated fixed points which build the hyperbolae C± re- 
ported in Table I. Given that their interpretation involves 
technicalities which do not add much to the current dis- 
cussion we submit the reader to [ll| for further details. 

Summarizing, in this subsection we have reinforce our 
results in pdT |. as we have shown that tracking solutions 
can exist in quintom cosmologies. We have identified 
that (provided they exist) these tracking solutions can 
be found in the region <fi — ±oo and tp = ±oo. 9 In 
addition, we have shown that the region 4> = ±oo is the 
only possible locus for the phantom cases. 



9 This feature was first noticed in JStJ but in the framework of 
single scalar fields. 



V. CONCLUSION 

The dynamical systems analysis is a very powerful tool 
when one wants to extract useful information about the 
dynamics of the cosmological evolution without knowl- 
edge of exact solutions and of specific initial conditions. 

In the present paper we have applied dynamical sys- 
tems tools to retrieve information about the dynamics 
of quintom models with arbitrary potentials (other than 
the exponential ones). Appropriate phase space variables 
has been chosen so as to allow for a study of regions in 
the (0,</?)-plane for both finite and infinite values of the 
scalar field variables cp, ip. 

The results of the present study indicate that existence 
of phantom attractors is not generic; it depends on the 
features of (the natural logarithm of) the potential. For 
quintom models, in the finite (</>, (^-region, there may ex- 
ist either de Sitter attractors associated with the saddle 
points of In V, or tracking attractors in the infinite ((/>, 
(^)-region. For the example studied in section III C (po- 
tential the de Sitter phase represents an attractor 
solution of the phase-space 

We want to underline that the present study is not 
particularly well suited for potentials which are products 
of sums of exponentials on each of the scalar fields. The 
reason is that in the latter cases dimensional reduction 
of the dynamical system is possible (and performing it 
is strong advisable). Therefore, our study should be re- 
garded as complementary to that of references d, [l(| HH • 
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APPENDIX A: DYNAMICAL SYSTEMS 
RECIPES 

The first step to obtain qualitative information about 
the solutions of a differential equation x' = f(x) is to 
study the flow of the differential equation in the vicin- 
ity of the critical points, i.e, to study their stability [281 ] . 
The essential idea is to linearize the differential equation 
at each fixed point and then use the Hartman-Grobman 
theorem which in colloquial terms can be stated as fol- 
lows: the flow of a non linear differential equation in the 
neighborhood of a non-hyperbolic critical point can be de- 
formed continuously into the flow of its linearization; i.e, 
the orbits of both systems are (locally) qualitatively the 
same. 

Since the flow of a linear differential equation can be 
obtained explicitly, we will present the analysis for linear 
systems. 

For the linear differential equation x' = Ax defined in 
K™ (or defined in a subset of the whole real space; or in 
general defined in a smooth manifold) we can determine 
the eigenvalues of the matrix A (complex in general and 
not necessarily different) and the associated eigenvectors 
which span three subspaces of R n (which depend on the 
nature of the associated eigenvalues): E s , E u y E c . The 
former subspaces contain disjoint orbits which form a 
partition of the state space, i.e. E s ® E u ® E c = W 1 . 
They are, respectively, the stable subspace (spanned by 
the eigenvectors whose associated eigenvalues has nega- 
tive real parts), the unstable subspace (spanned by the 
eigenvectors whose associated eigenvalues has positive 
real parts), and the central subspaces which are gen- 
erated by the eigenvectors whose associated eigenvalues 
have zero real parts. 

The stable and unstable subspaces are characterized, 
respectively, by the properties: 

x € E s => lim e Ar x = 0, 

r — >oc 

xeP^ lim e Ar x = 0. 

r — > — oo 

These describe the asymptotic behavior: all the initial 
states in the stable subspace are attracted by the critical 
point x = and all the initial states in the unstable 
subspace are repelled by x = 0. 

If the system x' = f(x) is non linear, we can use the 
Hartman-Grobman theorem. Using it we can define the 
manifolds £( s > u ' c ) (stable, unstable and center manifold 
respectively) at a fixed point, those manifolds are tan- 
gent to the corresponding subspaces E^ s ' u ' c ^ of the lin- 
earization at the fixed point (stable, unstable and center 
subspaces respectively). All the orbits in £ s converges 
asymptotically to the fixed point when the time variable 
evolves (t — > oo), whereas all the orbits at £ u converge 
asymptotically to the fixed point when r — > — oo. The 
manifold £ c contains all the orbits whose asymptotic be- 
havior can not be analyzed making use of the linear anal- 
ysis. 



APPENDIX B: APPLICATIONS 

Our dynamical systems were constructed making use 
of a maximum of five dynamical variables x l , i = 1 . . .5. 
The corresponding phase-space equations can be written 
symbolically by 

x*' = f(x\x 2 ,...,x 5 ), (Bl) 

where prime denotes the derivative with respect to the al- 
ternative time variable r which is chosen by convenience. 

The next step in the study of the evolution of our dy- 
namical system is to find its fixed (or critical) points 
(x 1 * ,x 2 * , . . . , x 5 *), which are given by the conditions 

/'•:.'• <•'": 0. (B2) 

The stability of the fixed points (x 1 * , x 2 * , . . . , x 5 *) is then 
analyzed by studying the linearized dynamical system ob- 
tained by expanding the evolution equations about those 
fixed points (as explained in many seminal references, 
c.g [30]). After that, one tries solutions in the form 
(x 1 , x 2 , . . . , x 5 ) = [c-x, C2, . . . , C5) e xt in the linear approx- 
imation, and finds that their characteristic exponent A 
and the constant vector (ci, ca, . . . , C5) must be respec- 
tively an eigenvalue and an eigenvector of the matrix 



dx 1 ' 


dx 1 ' 


dx 1 ' 


dx 1 


dx 2 


dx 5 


dx 21 


dx 2 ' 


dx 2 ' 


dx 1 


dx 2 


dx 5 


dx 5 ' 


dx 5 ' 


dx 5 ' 


dx 1 


dx 2 


dx 5 



(B3) 

The character of the fixed points depends on the val- 
ues of the characteristic exponents: if the real part of 
all characteristic exponents is negative, the fixed point 
is asymptotically stable, i.e., an attractor. On the other 
hand, it is enough to have (at least) one characteristic 
exponent with positive real part to make the fixed point 
asymptotically unstable (commonly called as source or 
repeller). The repeller is a saddle point if at least one 
of the other characteristic exponents has a negative real 
part, in which case, apart from the unstable manifold, 
there is a stable manifold containing the exceptional or- 
bits that converge to the fixed point. 

In addition, when one of the exponents is null the point 
is not hyperbolic and therefore structural stability can- 
not be guaranteed (the geometric form of the trajectories 
may change under small perturbations). Hence, the case 
in which the largest real part is precisely zero must be 
analyzed using other methods. In this case the linear 
analysis is inconclusive (the Hartman-Grobman theorem 
fails) . 
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